O 
O 



c3 



ANALYTICAL RAMIFICATIONS OF DERIVATIVES VALUATION: 
ASIAN OPTIONS AND SPECIAL FUNCTIONS 

Michael Schroder 

Lehrstuhl Mathematik III 
Seminargebdude A5, Universitdt Mannheim, D-68131 Mannheim 



Averaging problems are ubiquitous in Finance with the valuation of the so-called Asian 
options on arithmetic averages as their most conspicuous form. There is an abundance 
of numerical work on them, and their stochastic structure has been extensively studied 
by Yor and his school. However, the analytical structure of these problems is largely 
\ unstudied. Our philosophy now is that such valuation problems should be considered 

r -|r as an extension of the theory of special functions: they lead to new problems about 

new classes of special functions which should be studied in terms of and using of the 
methods of special functions and their theory. This is exemplified by deriving integral 
CN . representations for the Black-Scholes prices based on Yor's Laplace transform ansatz to 

their valuation. They are obtained by analytic Laplace inversion using complex analytic 
methods. The analysis ultimately rests on the gamma function which in this sense is 
found to be at the base of Asian options. The results improve on those of Yor and have 
served us a as starting point for deriving first time benchmark prices for these options 

1. Introduction: This paper studies analytical ramifcations of derivatives valuation 
which it considers as an extension of the theory of special functions. More precisely, the 
focus is on the analytical aspects of averaging problems. These are ubiquitous in Finance 
and make their appearance in as diverse streets cts interest rate derivatives, credit risk, and 
problems originating in life insurance contracts. We consider their most conspicuous form: 
the so-called Asian options on arithmetic-averages. These are widely traded options whose 
prices, in a way made precise in §4, are quoted using mathematical models. There is by 
now an abundance of numerical work on how to compute these prices with almost daily 
additions like [V]. However, even in the rather simplifying but still basic Black-Scholes 
model reviewed in §3, the conceptual understanding of their valuation appears to be far 
from simple and has lead early researchers to conclude that 'it is impossible to derive an 
explicit analytic expression for an Asian option price'. And as it can be seen from the 
discussion in [RS] for instance, coming to grips with valuing Asian options has intrigued 
financial theorists for over a decade by now. 



O 

a 

> 

An explanation for this has in a sense been found by Yor and his collaborators whose work 
\ has defined the state of the art on this problem in particular. Establishing connections with 
Bessel and confluent hypergeometric functions they have studied the stochastic structure 
of Asian option valuation with [Y01] as the lastest addition. With regrad to closed form 
valuation formulas for Asian option prices, Yor defined in [Y] the standard with his triple 
integral. And as a second line of thought, the Laplace transform approach developed 
in [GY] was a celebrated advance. Its results, and its actual relation to Asian option 
valuation, is reviewed in §8 and furnish the starting point of our research. 

The main mathematical result of this paper is a second valuation integral for the Black- 
Scholes price of Asian options which improves on Yor's. As described in §6 it expresses this 
price as a single integral over now, roughly speaking, the product of two well-studied higher 
transcendental functions. Both are given as integrals and built up using so-called Hermite 
functions, reviewed in §2. These last functions generalize the familiar complementary error 
function and come from boundary-value problems in potential theory for domains whose 
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surface is an infinite parabolic cylinder. It seems to be the first time that such parabolic 
cylinder functions are explicitly identified to characterize solutions to problems of Finance. 
Moreover, our results improve on those difficulties with numbers of gigantic sizes which 
render Yor's triple integral so uncomputable. This we discuss in §7. And, in fact, our 
formula improves on these difficulties to such an extent that it has served in [SH] and 
[ST] as the starting point for deriving first time explicit benchmark pricing formulas for 
Asian options. 

Mathematically, our valution formula is obtained by analytic Laplace inversion of certain of 
Yor's Laplace transforms. And this can be regarded as the main mathematical contribution 
of the paper. The details of the proof are described in Part III of this paper. Its methods 
are complex analytic. It uses the complex inversion formula for the Laplace transform and 
is ultimately based on the classical Hankel formulas for the gamma function. However, 
inversion is achieved under restrictions only that are of stochastics origin. These are then 
lifted in a second step using analytic continuaton. This seems to be a new feature which 
illustrates how complex analytic methods can be used to extend the validity of stochastics 
results in situations where stochastic methods alone seem to lead to intractable problems. 

And all this illustrates part of our guiding philosophy that such valuation problems seem 
best considered as extensions of the theory of special functions. They give rise to new 
problems in and new classes of special functions which should be studied in terms of and 
using the methods of special functions and their theory. 
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2. Preliminaries on Hermite functions: This section collects relevant facts about 
Hermite functions from [L, §§10.2ff]. These functions are given in terms of the Kummer 
confluent hypergeo metric function $, and the Hermite function of degree any complex 
number \i is the function on the complex plane defined by 



Hermite functions so specialize to the /x-th Hermite polynomials if /i is any non- 
negative integer. Moreover, they are holomorphic on the complex plane as functions of 
both their variable z and their degree \i. If \i is not a non-negative integer, they have the 
absolutely and compactly convergent series 
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Preliminaries and recollections 
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which also describes their behaviour for small arguments. If the real part Re (/z) of [i is 
negative, the asymptotics of for large arguments is described by the expansion 



valid for any complex z with positive real part. Herein recall (a)k as the Pochhammer 
symbol given by (o)o = 1 and (a)k+i = (a+k)(a)k for k any non-negative integer. 

To discuss further connections of Hermite functions with other classes of special functions, 
if the real part of fi is negative, we have the integral representation: 

r>00 



i r 



for any complex number z. Hermite functions thus also generalize the complementary error 
function Erfc which for any complex number z is given by: 



9 r°° 

Erfc (z) = e~ z2 -= / e~ u2 - 2zu du ; 
V 71 " Jo 



in fact, (2/y / 7r)i7_i(^) = exp (,2 2 )Erfc (z). Hermite functions of degree any complex 
number \i are moreover connected with the parabolic cylinder functions by D^(z) = 
2 _ m/ 2 ex p (—z 2 /4)H^(2 1 / 2 z) for any complex z, and with the second Kummer confluent 
hypergeometric function \1> by H^z) = 2 M \I/(- / u/2, 1/2; z 2 ) if Re (z) > 0. 



3. Black— Scholes modelling: The analytical problems to be discussed originate from 
the so-called risk-neutral approach to the valuation of contingent claims. Textbook treat- 
ments for this and other notions developed for the analysis of financial markets and in- 
struments are in [DSM], [MR], [KSb, Chapters 1-4] for instance. This analysis is based 
on modellings of security markets, and this section aims to give some background on the 
most fundamental of these, the Black-Scholes modelling of security markets. 

In fact, we need only that particular case of the Black-Scholes model where there are only 
two securities., and the understanding is is that these are traded on markets where their 
prices are determined by equating demand and supply. First, there is a riskless security, a 
bond, whose price (3 grows at the continuously compounding positive interest rate r, i.e., 
for which we have we have f3 t = exp (rt) at any time t G [0, oo). Then, there is a risky 
security, and the fundamental idea is that all uncertainties affecting its price £ yield a 
certain probability space. In fact, consider for this a complete probability space equipped 
with the standard filtration of a standard Brownian motion on the time set [0, oo). Giving 
expression to the fact that 5 comes as an equilibrium price, we have the risk neutral 
measure Q on this filtered space, a probability measure equivalent to the given one. And 
with B any standard Q-Brownian motion, the exact modelling then is that S is the strong 
solution of the following stochastic differential equation: 



dS t = vj- Sfdt + a- Sf dB t , t e [0, oo) , 
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or equivalently using ltd calculus, 

St = s e im -^ 2)t+aBt te [0,oo) 



The positive constant a is the volatility of S. The specific form of the otherwise arbitrary 
constant w depends on the nature of the security modelled which could be a stock, a 
currency, a commodity etc. For example, if S is a stock paying a dividend at the continuous 
rate S, we have w = r — 5. 

4. Asian options and their equilibrium pricing: In the Black Scholes framework 
of §3, fix any time to and consider the accumulation process J given for any time t by: 

J it) = S u du 7 

The arithmetic-average Asian option written at time to with maturity T and strike price 
K is then the stochastic process on the closed time interval from to to T paying 



T-t ) { T-t 

at time T and paying nothing at all other times. As such it is a contingent claim on the 
time interval from t to T with payoff (J(T)/(T-t ) - K) + . 

It is one of the fundamental insights that in the equilibrium framework of the Black-Scholes 
model any such contingent claim on a risky security has an equilibrium price too which is 
equal to the expectation of its payoff with respect to the risk neutral measure conditional 
on today's information, see [MR, Corollary 5.1.1] for instance. Applying this arbitrage 
pricing principle, the price Ct of the Asian option at any time t between to and T so is 
given as the Q-expectation conditional on the information available at time t: 



C, = e-^-OE" 



AT) K 



T-t 



However, following [GY, §3.2], do not focus on this price, but normalize the valuation 
problem as follows. On factoring out the reciprocal of the lenght T— to of the time period, 
split the integral J(T) into two integrals, one of which is deterministic by time t and the 
other of which is random. Couple the deterministic integral with the new strike. For the 
random integral, restart the Brownian motion driving the underlying at time t, and then 
using the scaling property of Brownian motion, change time to normalize its coefficient in 
the new time scale to two. The precise result is the factorization: 



T-t ' a 

which reduces the general valuation problem to computing 

cM(h,q) = E°[(AM- q )+], 
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the normalized time-t price of the Asian option. To explain the concepts, is Yor's 
accumulation process 



A W) = f e 2(B w +urv) dW) 

Jo 



and the normalized parameters are as follows: 



where 



z/=^-l, h=^(T-t), q = kh+q*, 
k=f^ q* = q*(t) = ^r( K - (t-*o) ~ J* S u du^J . 



To interpret these quantities, v is the normalized adjusted interest rate, h is the normalized 
time to maturity, which is non-negative, and q is the normalized strike price. 



Part II Statement of results 

5. A reduction of the valuation problem: Computing the normalized time-t price 
of the Asian option reduces to the case where the normalized strike price q is positive. 
Indeed, if q is non-positive, Asian options loose their option feature, and their normalized 
time-t price is given by 

C<"> (h, q) := E Q [(A™ - q) + ] = E Q [A™] - q . 
On applying Fubini's theorem, this last expectation is computed as follows: 

r < m e 2h{v+1) - 1 
^ PJ" 2(z/+l) ' 

for v any real number, and being equal to h if v = — 1 in particular. 

6. Statement of the main results: The main result of this paper is a closed form 
solution for the normalized price (h, q) of the Asian option in the generic case where 
q is positive. It expresses this function as the sum of integral representations. These 
are obtained by integrating the product of Hermite functions discussed in §2, with 
functions derived from weighted complementary error functions. Proved later in Part III, 
the precise result is as follows: 

Theorem: If q is positive, the normalized price C^ u \h,q) of the Asian option is given 
by the following difference 

C^(h,q) = ce 2h ^S u+2 -cS u 
where the S% are three-term sums 

S$ = C tri g,e 5 ,c(P?) + Chyp,0 € ,£(p$) + Chyp,0 4 ,-e(P€) 

whose single summands are integrals that depend on 0% G [§ ,7r] and > 0, but which as 
a whole are independent of these parameters. 
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To explain the concepts, c is given by 

I>+4).(2g)^+ 2 ) 
c = c (v, q) = j — 

recalling z/ = 2vj/a 2 — 1. With 6* any real in the closed interval [-|, 7r] and p any non-negative 
real, the trigonometric terms Ctrig,6>,£(p) with £ equal to v or z/+2 are the integrals 

GWP) = ^Re ^-(,+4) ( - ) W (P+^)) # 

over the real parts of products of Hermite functions times certain functions E b (h), and the 
hyperbolic terms Ch yp ,e,^(p) with £ equal to ±z/ or ±(z/ + 2) are the integrals 

C hypA dP) = (g-(, + 4)(- C ° Sh Jg^ ) )^(^)(y+^))^ 

over the imaginary parts of such products. Herein E^(h) are the weighted complementary 
error functions for any complex w given by 

E^h)(w) = e^Erfc + ^v^) . 

Remark: If p equals zero, the trigonometric terms specialize to 

Ctri g ,^(0) = 2^H_ {v+4) (" £ ^) cos(^)#. 



Notice that for iH-4 any non-positive integer, the Hermite functions of the formula specialize 
to the corresponding Hermite polynomials. As explained in [SH, §7] and [ST, §4 Corollary], 
the normalized price is then given by a finite number of weighted complementary error 
functions and their derivatives. This effect is also observed in [Du] where Asia densities 
are expressed as certain integrals against Hermite functions. 

Finally deriving benchmarks for the normalized prices of Asian options seems to be one of 
the main practical application of formulas as that of the Theorem. To get an impression, 
Asian option on stocks are usually written with maturities of one year and issued at par: 
K = So taking to = 0. Current annual volatilities for stocks range between 20 and 50 
percent. Taking an annual interest rate of nine percent, we then have from [ST, §12] 



a 


maximal error 


CW(/i, q) 


20% 


4.9727 x 10" 16 


. 00074155998788343 


30% 


4.9687 x 10" 16 


0.00217354504625037 


40% 


4.9157 x lO" 16 


. 00478100328341654 


50% 


4.9461 x 10" 16 


0.00890942045213227 



Table 1. Normalized prices C (u) {h,q) of the Asian option for T=l. 
as benchmark normalized prices for the Asian option sharpening the results of [RS] . 
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7. Yor's triple integral and a comparison: The closed form results of §6 have all 
been preceded by Yor's integral representation [Y, (6.e), p. 528]. For the normalized time-t 
price (h, q) the latter implies the following representation as triple integral: 



CH(/i, q) 



setting 



x 



2 



--q 



i>xy{h) dydx 



1 



7t_ 

,2h 



and where for any positive real number a, the function ip a is given for any h > by 

-a cosh(w) 



Mh) = f 
Jo 



ml 

" 2h . 



sinh(w) • sin 



—w) dw . 
h J 



While Yor's formula seems to require v to be bigger than at least minus one, it is valid for 
all reals v. This is proved in [SL, §6]. Our formula, in contrast, is given as a sum of single 
integrals whose integrands have a structural interpretation as products of two functions. 
It identifies the higher transcendental functions occuring as factors in these products, and 
shows how they are given by or built up from Hermite functions. On a technical level, these 
differences can be regarded as consequences of the different mathematical approaches for 
proving the valuation formula. In fact, Yor's originates from a direct attack on the Asia 
density whereas our's is eventually based on the indirect enveloping construction of [CS], 
and thus on non-Asian options. 

However, there are not only structural differences between the two results. One purpose 
of explict pricing formulas is to provide means for actually computing option prices in real 
life. And here Yor's formula has a number of structural difficulties. For instance in the 
setting of the §6 example we compute the factors c Vj h as follows: 



Cv,h 


20% 


a = 30% 


a = 40% 


a = 50% 


T = 1 year 


4.379 x 10 214 


2 . 647 x 10 95 


4 . 266 x 10 53 


1 . 753 x 10 34 


T = 6 months 


1.321 x 10 429 


6.717 x 10 19U 


2 . 288 x 10 1U7 


4 . 846 x 10 s8 



Table 2. 



C v ,h 



as function of T and u. 



Yor's formula thus expresses the normalized price of Asian options, which is not too big, 
as the product of a big number times a triple integral. The latter so has to be small. And 
it must be computed with very high accuracies to get reasonably accurate results. The 
formulas of §? Theorem seem to be better behaved in this respect. For instance, in the 
case T = 1 and a = 30% we have from [SH, §11] as values for the hyperbolic terms 







V 


-1534170.07497 


—v 


1534169.98118 


z/+2 


-4198799.82516 


-(I/+2) 


4198799 . 74048 



Table 3. Values of c-C hyPt *^(0) for T=l and CT =30%. 

where c = 4.072 x 10" 12 . And using any of the formulas of [ST, Part II] we get 0.00104 
as total contribution of the trigonometric terms. 



8 



Michael Schroder 



Part III Proof of the valuation formula 

8. Reduction to the Laplace inversion of certain Laplace transforms: The 

Laplace transform has been forged by Yor and his school into a very effective tool for 
transferring stochastics problems into analysis, and [Y01] is the last addition to his work 
on averaging problems. It seems to be characteristic for these problems that it is Bessel 
function which enter at the Laplace transform level. While this is briefly reviewed too, 
this section's aim is more precisely to reduce computing the normalized time-t price 



C (u) := E Q 



(A^-(kh+ q *)) + 



of the Asian option introduced in §4 to the following Laplace inversion problem. 

Proposition: If the positive real number h is such that q(h) = kh + q* is positive, the 
normalized price of the Asian option at h is given by: 

CM = i? -i(> Gr (<?(/*), *))(/*) 

as the Laplace inverse of FGY(q(h), z) at h. 

Here we use the concepts of §4. To explain the functions, for any real a > define: 

D v (a,z) 



F GY (a,z) 



z- (z-2(v + l)) ' 



for any complex number z with positive real part bigger than 2(v+l), where on choosing 
the principal branch of the logarithm: 

D v (a,z)=^- [ e-^-x^-I^-^(-) dx . 



a Jo ^va 

Here 1^ is, for any complex number //, the modified Bessel function of order \i more fully 
discussed in §7 below or in [L, Chapter 5], for instance. 

The basic idea originating with Yor for proving such results is to make the time h a variable 
and compute the Laplace transform of the functions so obtained. In the case of CM this 
is made difficult by the time variable entering not only via the stochastic process but also 
via the normalized strike price. Following [C] , as a first step in proving the Proposition so 
consider for any real number a the non-Asian option price functions foY,a on the positive 
real line that send any x > to 

f GY , a (x) = E Q [(A^-a) + }. 

Taken individually they cannot be used to value the original Asian option. However, as a 
whole they allow one to recover CM. Indeed, if q = kh+q* is positive, CM is obtained by 
choosing the function fcY,kh+q* and evaluating it at h. Equivalently, using the injectivity 
of the Laplace transform on continuous functions, CM so is equal to the Laplace inverse 
at h of the Laplace transform of fGY,kh+q*- Thus we are reduced to identify any FcY,a as 
the Laplace transform J^(fGY,a) of foY,a- 
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For explaining how Bessel functions enter into the picture, we restrict to the case v > 
where S'ijGY.a) has been computed in [GY, §3]. This will illustrate the essential stochastic 
ideas, and is also sufficient for our purposes. A proof for arbitrary v using complex analytic 
tools in particular can be found in [CS]. The punchline of Yor's computations is as follows. 
The basic idea is to make time stochastic using the Lamperti identity 



for v > proved in [Y92a, §2]. Herein w- u > is the Bessel process of index v and with 
_R(^)(0) = 1, see [RY, XI §1]. Then, has the double role of both underlying and 
stochastic clock. Using the fact that a is independent of time, thus transcribe the condition 
on the underlying to be bigger than a as the first passage time 

r Vi a = mf{w\A%> > a} 
for the stochastic clock. This is the key idea for obtaining the representation 



f GY ,aH = E Q 



2(iH-1)[w-t„, ] + 



1 



for all w > 0. At first sight this may have complicated the problem. However, it is just 
what fits the Laplace transform L of foY,a- 



L(z) 



f 

Jo 



> E Q 



e 2(i/+l)[tu-T„,„] + _ 1 

2(z/ + l) 



dw 



For computing this integral to show L = Fcy,a it seems best to follow a communication 
of Yor's, switch to measurable functions and applying Tonelli's theorem interchange the 
Laplace integral with the expectation E®. In this way we get the identity 



L{z) = 



1 



z(z-2(u+l)) 



E Q 



of measurable functions for any complex number z with Re (z) > 2(u+l). To identify the 
expectation in the numerator as D v a write it out as 



x 



. E Q 



X 



■Pu,a{x) dx 



where p v ^ a is the Bessel semigroup of index v starting at 1 at time a. Explicit expressions 
for them are classically known for v > —1, see [Y80, (4.3), p. 78] or [GY, Proposition 
2.2], while the case v < — 1 has been addressed only recently in [YG6, §3]. Under the 
hypothesis v > 0, Yor has in [Y80, Theoreme 4.7, p. 80] (see also [GY, Lemma 2.1 and 
Proposition 2.6]) explicitly computed the conditional expectation factor of this integrand 
while for v < there are no results. The upshot is that these results are available both 
only if v > 0, and this is where the restriction on v comes from. On substitution we get 



D„(a, z) 



l 

la 



a 



x 2 
"Ta 



x" +3 I V2:Hv2 



X\ 

a J 



dx , 
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as required. To complete the Tonelli argument proposed by Yor and to complete the proof, 
we have to establish the finiteness of this integral for any fixed complex number z with 
Re (z) > 2(z/+l). There is a further technical point to be taken care of. Choose the 
principal branch of the logarithm to define the square root on the complex plane with the 
non-positive real line deleted, whence \x = (2^+z/ 2 ) 1 / 2 has a positive real part. Integrability 
now follows using the asymptotic behaviour of the Bessel function 1^ near the origin and 
towards infinity, and we have shown L = Foy,a if ^ > as required. 

9. Preliminaries on integral representations of gamma and Bessel functions: 

With hindsight, it is the gamma function which is at the base of valuing the Asian option. 
However, the gamma function does not enter via Euler's integral recalled to be given by 

poo 

r(s) = / e~ x x s ~ 1 dx, 
Jo 

for any complex s with Re (s) > 0. The problem here is that such a formula would be 
required for s with non-positive real parts where the resulting integrand is not integrable 
near the orgin. Being able to perform such integrations from minus infinity to plus infinity 
avoiding the origin needs higher dimensions, and recall Hankel found natural candidates 
for this some 130 years after Euler. His Hankel contours Cb,r with parameters 9 any angle 
in (f, 7r) and R > any real are the following contours of integration in the complex plane. 
One comes in from infinity on the ray {yexp(— i9) : y > R} until the point i?exp(z6*), 




Figure 1. The Hankel contour Cq^r. 



then passes counterclockwise around zero with distance R until the point _Rexp(i6>), and 
leaves on the ray {yexp(i9) : y > R} to infinity. And using integration over any of these 
contours, we then have HankeVs formula: 



1 



T(s) 2m 

valid now for any complex number s, see for instance [D, pp.225f]. 

The occurrence of modified Bessel functions in §6 Proposition now furnishes the point of 
attack. From [L, §5.7] they are given for any complex number fx by the following series: 



WJ Z^ m=0 mirGu + m + 1) ' \2J 



for any complex z in C\R<o- The idea is to substitute for any of the reciprocal gamma 
factors of the single terms of this seris substitute the corresponding Hankel formulas and 
interchange the order of summation and contour integration. With the details in [WW, 
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17-231, p. 362], for instance, the result is the following Hankel-type integral representations 
of modified Bessel functions on the right half-plane 

Lemma: For any modified Bessel function I M one has: 

Z7rz JlogC 

for any complex number z with positive real part and any Hankel contour C . 

Herein the principal branch of the logarithm on the complex plane with the non-positive 
real axis deleted has been chosen. The contour logC^ thus has the following shape: 
Coming in from plus infinity, move on the parallel through the point —id to the real axis 

log R+i9 





1 s> 


logfl 


1 >■ 

1 < 



log R-id 

Figure 2. The contour log Cq,r- 

to the point \ogR — id. From this point move up to the point logR + id on a parallel to 
the imaginary axis. Finally exit from \ogR + id to plus infinity on the parallel through id 
to the real axis. Notice logC^ C {z\Re (z) > 0} if and only if the radius of the circle in 
C e , R is bigger than 1. 

10. First steps of the Laplace inversion: This section resumes the proof of §6 Theo- 
rem by addressing the Laplace inversion problem of §8 Proposition. Define for any positive 
real a and for b equal to v or u + 2 functions on the positive real line by sending any h to 



Z a b (h) = [°°— [ e -* 2+ *W COsh (™V+3,W 6 ~ WVT \ dwdx. 



1 f - x 2 +x ^ cos h(w) e~ w ^ 

' log C s ,r 

where log Cq^r is any logarithmicalized Hankel contour as in §7, and let c\ denote 

-Wh (2a) 
ci = e • -i — '- . 

(z/ + l) • e 2 ^ 

Then §8 Proposition reduces to compute the functions as follows 

Lemma: If v is non-negative, the inverse Laplace transform at any positive real number 
h of any function -F^y (a, z) is given by: 

^- 1 (F Gy (a,-))(/ i ) = c x {r, a v+2 {h) -Kih)) , 

where any of the functions is independent of the Hankel contour chosen and has an 
absolutely integrable integrand. 
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The proof of the Lemma is based on the complex inversion formula for the Laplace trans- 
form, see [B, Chapter 7, §5]. 



To prove the Lemma, first recall that for any positive real number a we have: 

F GY (a, z) = — o: , 

z ■ [z - 2(u + l)j 

where 

D u {a,z)=—- \ e~ia .x^.l^-^l \ dx, 

for any complex number z with real part bigger than 2(H-1) . Then let z' be any positive real 
number such that the line {z|Re (z) = z' } is contained in the half-plane where Fqy{cl, z) 
is a holomorphic function. Suppose for a moment proved that this function is rapidly 
decreasing, and apply the complex inversion formula to it. Writing out D l/ (a,z), the 
inverse Laplace transform of Fqy{cl, z) is the function on the positive real line given for 
any positive real number h by: 

x 2 

1 rz' +ioo i-oo ~2a . T .^" 1 " 3 . T ( X\ 



e hz , J^ {a) dxdz. 

2™Jz'-ioo Jo z(z -2(u+l)) 



For an equivalent expression, change variables rj = 2z+v 2 , put zq = 2z' +v 2 , and substitute 
the Hankel-type integral representation of §9 Lemma for the modified Bessel function. 
With c[ = 4(u + 1) ■ (2a)~^ +4 )/ 2 • c\, the following integral is then to be computed: 



c\ ( | / [ [ e 2a x +3 e a • e 2 ; -rrrdwdxdz 

1 \2ixi) J zo _ loo J J logC (z - u 2 )(z - (v+2) 2 ) 



We claim that for R sufficiently big the absolute value of the integrand of this integral is 
exponentially decreasing to zero with the absolute values of z, x, or w going to infinity. 

Granting this result, the above triple integral then gives the desired Laplace inverse, and, 
using Fubini's theorem, the order of its integrals can be interchanged. Take the integral 
for Laplace inversion, i.e., the integral over the line {^|Re(^) = zq}, as inner integral. 
Change variables x = (2a) 1 / 2 t. The Lemma follows on decomposing the denominator of 
the integrand in partial fractions. 

One is thus reduced to prove the last claim about the asymptotic behaviour. For the 
calculations recall | exp (£)| = exp (Re (£)), for any complex number £. For the asymptotic 
behaviour in the absolute value of w, reduce to elements w = x ± i9 in \ogCe,R- For these 
Re(cosh(w)) = — cosh(x)| cos(#)|. The absolute value of the hyperbolic cosine factor of 
the numerator equals exp ((x/a)Re (cosh(w)), whence the required asymptotic behaviour 
in the absolute value of w. 
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For the asymptotic behaviour in z, notice | exp (— wz 1 ^ 2 )] = exp (— Re (wz 1 ^ 2 )) and recall 
z 1 ! 2 = exp ((1/2) • (log \z\ + i arg(z))). The argument of z converges to n/2 with \z\ going 
to infinity. Now Ke(wz 1 / 2 ) equals \z\ 1 ^ 2 times Re (w) cos(arg(z)/2) — Im(w) sin(arg(^)/2). 
Herein the cosine is positive and bigger than cos(7r/4). Thus Re (wz 1 / 2 ) is positive if the 
real part of any w is big enough. Hence choose R big enough to have the desired asymptotic 
behaviour in the absolute value of z. 

The asymptotic behaviour in x is determined by exp (— x 2 /2a), thus completing the proof. 

11. Computation of certain Laplace transforms: This section computes the in- 
verse Laplace transforms identified in §10 Lemma thereby illustrating typical techniques 
for working with Laplace transforms. 



For any a and (3 in C, consider the functions on the positive real line given by: 

9<.A*) =4" ( e+ ° Erfc (^r + pVr ) + e "°'" Erfc (^T " /3VF )) ' 

for any positive real number t. Then one has the following two results: 
Lemma: If the real parts of a and a 2 are positive, one has: 

for any complex number z in C \ R<o with real part bigger than |Re (/3)|. 
Corollary: If the real parts of a and a 2 are positive, one has: 

2* (g a ,p) (z) = - , 
for any complex number z in C \ R<o with real part bigger than Re (f3 2 ). 

The Corollary follows from the Lemma upon decomposing the denominator in partial 
fractions and using the linearity of the Laplace transform. 

We use the following two results proved mutatis mutandis in [D, Beispiel 8, p.50f]: 

2 7 {iP a ){z) = e- a ^ where ^ a {t) = — = • t~* ■ e , 



1 — 



(I 



2 



&(Xa)(z) = — where x<x(t) = -7= -t 2 • e 

\/Z V 7T 



for any complex number z in C \ R<o and any positive real number t. 
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Subtracting Xa from f a ,/3, the proof of the Lemma reduces to show the identity: 

S (-„ • E* (W + ^)) W = -/J^j^. 

Multiplying any nice function / with exp (a-) induces a shift by —a in its Laplace transform: 
(exp (at)f(t))(z) — J^(f)(z — a). Using this with a = (3 2 one is further reduced to 
calculating the Laplace transform of the above complementary error function factor only. 

Since Re (a) is positive, the real part of flu 1 / 2 + (a/2)u~ 1 ^ 2 goes to plus infinity with u 
converging from the right to zero. The Fundamental Theorem of Calculus thus gives: 

Using the transform-of-an-integral property 3? (J* f (u) du) (w) = S 9 (f)(w) the La- 
place transform of this complementary error function at w = z — (3 2 is given by: 

2/3 ^ ( 1 i -(PV^+7^r=Y\, n o, 

t ~2 . e v v 2V ^j (>-/3 2 ) 



2(z-p) \0F 



2 

-I 2- 



Using the Laplace transform of Xa, the first Laplace transform of this sum equals 



Using the Laplace transform of ip a , the second Laplace transform of this sum equals 



e~ af3 e~ a ^ . 



The above sum hence equals exp (—a/3) exp (— a^fz) /(z 1 / 2 (z 1 ^ 2 + 13)). The identity to be 
proved follows upon substituting this last expression. This completes the proof. 



12. Two intermediate results: In this section the Laplace inversion of §10 is resumed 
concentrating on the two summands of the expression derived in §10 Lemma. If v is 
non-negative, recall them as the functions on the positive real line sending any h to 

where a > and b are any reals, and logCg^ is any logarithmicalized Hankel contour as 
discussed in §9. For these integrals we have the following key result: 
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Proposition: If moreover R > 1, we have 
Et(h) = il>+4)e^ • JL J (E h {h) + E_ h {h)) (w) • H_ {v+A) (-^=$) dw 

where the respective integrands on both sides are absolutely integrable, and the integral on 
the right hand side does not depend on the Hankel contour chosen. 

Notice that the T, b are complex-valued integrals, and we so have the decomposition 

E a b (h) = Re(E a b )(h) +z!m(E a b )(h) 

into real and imaginary parts. Our second result then expresses the contour integral 
representation of the Proposition in terms of these real-valued integrals as follows 

Corollary: Under the additional assumption p = logR > 0, we have 

1 i u z* 00 

Re(£g)(*) = — J e"*V+ 3 J Re {F^ b (h))(p + icfy d<pdx 
Im(E£)(/*) = — e 5 ^ e-^x v+3 J 1m(F£ tb {h)){y+i6)dydx. 

To explain the notation in these two results, the Hermite functions of degree p are 
discussed in §2, and the two other functions that occur are given by: 

E b (h)(w) = e wb • Erfc ( + ^-V2h] , 

\V2h 2 J 

K,b( h ) H = e ^ C ° Sh(W) ( E b(h) (w) + E- b (h) («;)) , 
for any complex number w. 

Proof of the Proposition: For the proof of the Proposition choose a Hankel contour 
Ce,R with R so big that for any element w in logC^jj also the real parts of both w and w 2 
are positive. At any w in logCe^ substitute in .y b (h) the inverse Laplace transforms of 
§11 Corollary. Interchange the order of integration using the absolute integrability of the 
integrand. This gives the expression of the Lemma for ^f b {K). The integrand herein is a 
holomorphic function on C \ R<o- Using the Cauchy Theorem, the value of the integral 
is independent of the Hankel contour Cq^r chosen as long as log R is non-negative. This 
completes the proof. 

Proof of the Corollary: Given the Proposition, the proof of the Corollary is an exercise 
in path integration. Put F x = F£ b (h). Changing the order of integration there, 



oo 

2 



2™ J\o g Ce, R 



e~ x x u+6 — / F x {w)dwdx 
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is to be calculated. Concentrate on its inner integral and consider the following subpath 

P = P c + P ao with P c = P c ,+ - P c ,_ 

of the path logCg,/?. Here the path P c - starts from logR and moves parallel to the 
imaginary axis to the point log R — i9. The path P c ,+ starts from logR and moves parallel 
to the imaginary axis to the point log R + iO. The path P^ moves from log R + iO parallel 
to the real axis to +00. For any x > 0, the inner integral thus breaks up as follows: 

^- / (F x (w) - F x (w)^j dw + ^- [ F x (w) dw - ^- I F x (w) dw , 

" Poo " Pc , -\- " Pc , — 

with w the complex conjugate of w. Using the series expansion of the exponential and 
the complementary error functions, F x is compatible with complex conjugation, i.e., F x 
evaluated at the complex conjugate of any complex number w is the complex conjugate of 
F x at w: 

F x (w) = F x (w). 

Write the elements of P^ as w = y + id with y > log R, notice dw = dy and change vari- 
ables. Using the compatibility of F x with complex conjugation, one obtains the improper 
integral from logR to infinity of the imaginary parts of F x (y + i6). 

On the circle part P c of P one has to be a bit more careful about the volume forms. The 
elements of P c ,+ are parametrized by w = logR + i(p with <p in [0,0], whereas those of 
P c ,- are parametrized by logR — i(p with in [0, 9]. Changing variables accordingly and 
integrating from zero to 9, the induced volume form in the P c ,+ -integral is dw = i d(p, 
whereas that on the Pc, --part is dw = —i d(j). Using the compatibility of F x with complex 
conjugation, this completes the proof of the Corollary. 

13. Explicit calculations: This section explicitly computes the integrals of §12 Corol- 
lary in the case where the parameters p there are equal to zero. Recalling for this the 
function given on the positive real line by 

where a > and b are any reals and C = Cg t i is any Hankel contour as discussed in §9, 
the precise result to be proved for establishing §6 Remark for v non-negative is the 

Lemma: If v is non-negative, we have 

Re (EJ)(/i) = 2^ei> 2 " £ H_ ( „ +i) (-^) co sW )#. 
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The Lemma is proved by computing the integrals over F x := F® b of §12 Corollary upon 
choosing R = 1 there. For completeness sake start with the imaginary part integrals: 

POO 

/ Im (F x )(y + i9)dy. 
Jo 

Abbreviate j3± = y/V2h± (b/2)V2h, and write the complementary error functions oc- 
curing in F x (y + i9) as improper integrals starting from zero as in §2. Multiply out the 
expressions of the exponents to obtain: 

e ±( ^ )b Erfcf/3 ±+ z-L) 

Jo 

Since j3± is a real number and u can be taken as real numbers, the real respectively 
imaginary parts in F x are determined by the real respectively imaginary parts of the 
exponentials in the integral, whence in particular 



Be(E b (h))(y+iO) = 
Im (E b (h))(y+iO) = 




For proving the Lemma, 



/ Re^)^)^ 
Jo 



is calculated. In contrast to the above argument, in the complementary error functions 
occuring in F x (i<p) now the following paths of integration are used: Abbreviating (3± = 
±(b/2)y/2h, first move from (3± + i(fi/\/2h to f3±, then continue from (3± to plus infinity 
along the real line. Since cosh(z0) = cos(0), the above integral then equals: 

£ e K W- cos ^|(Erfc (J3+) + Erfc (/?_)) • cos (06) + -^($+ + #, 

upon abbreviating for any angle (f>: 



$ ± (0) = _Re (e ±i4 * J ± ^ e ~ u2 du^ 



In the first of these two last integrals notice Erfc (/?+) + Erfc (/?_) = 2 since (3- is minus 
/?+. To calculate $± change variables u = j3± + iw/V2h in the integral. Write the factor 
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i/y2h that is picked up as exp (z(7r/2))/V2/i . Multiplying out the expression obtained in 
the exponent, it follows 



J 8+ v2/i Jo 



>8± 

The real part of this expression is determined by the real part of the exponential functions 
in the integral. Abbreviating x u = b(9 — u), the values of the cosine at tt/2±x u thus appear 
as factors. A shift by tv/2 turns a cosine into a sine as follows: cos(7r/2 ± x u ) = ^fsm(x u ). 
Thus is minus $ + , and the expression of the Lemma results. 

14. First part of the proof of the valuation formula: The proof of the valuation 
formula of §6 Theorem is in two steps. As a first step, §6 Theorem is in this section 
established for v non-negative. As a second step, this equality is extended in the next 
section to any complex number v using analytic continuation. 

Thus let v be non-negative. Recalling §10 Lemma, the valuation formula of §6 Theorem 
is obtained by subtracting Ti%{h) from T,l +2 (h) and multiplying this difference with the 
constant c\. Substitute the expressions computed in §12 Proposition. Defining Sb by 

and cancelling the factors exp(/i^ 2 /2) thus completes the proof of the Theorem's repre- 
sentation for the normalized time-t price of the Asia option. 

15. Second part of the proof of the valuation formula: This second part of the 
proof of §6 Theorem extends its validity from v non-negative, as established in the previous 
section, to v any complex number. 

This reduces to show the following two results. First, the normalized price is an entire 
function in v. Second, the right hand side of §6 Theorem is a meromorphic function on 
the complex plane. Indeed, these two functions agree on non-negative real numbers v. 
Using the identity theorem they so agree on the complex plane as meromorphic functions. 
With one of them entire, the other one is then entire, too. In fact, for establishing the 
second result it is sufficient to divide off the function c with its at most simple poles in the 
negative integers and prove the following result: The functions Sb of the right hand side 
of §6 Theorem are meromorphic with at most simple poles off the negative integers. 

For proving the first of these two results recall the normalized price C^(h, q) as given by 

C^(h,q) = E[f^)(h)J 

where A^ v \h) = f Q h exp ((2(B u +uu)) du, and we have f(x) = (x—q) + , for any real x. The 
proof of this normalized price being entire in v is then further reduced to show 



E f(A h )e 



an entire function in v where A is the process A(°). Indeed, this follows using the Girsanov 
identity E[f(AM(h))] = exp (-u 2 h/2) E[f(A(h)) exp (uB(u))] of [Y, (l.c), p.510]. 
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For proving the above function entire develop the factor exp (uBh) of f(Ah) exp (uBh) into 
its exponential series. Suppose computing the expectation of the resulting series term by 
term is justified for any complex number v. For any complex number v, this then gives a 
power series in v and thus explicitly shows E[f(A h )exp (uBh)] holomorphic at v. 

Interchanging the order of integration and summation is justified using Lebesgue Domi- 
nated Convergence if the following is true. The expectations of the absolute value of any 
single term of the series for f(Ah) exp (uBh) exist and the series so obtained converges. 
Using the Cauchy-Schwarz inequality this is implied by the series: 



being convergent for any complex number v. Implicit herein is that j{Ah) is square 
integrable. This is implied by Ah being square integrable. Any n-th moment in particular 
of A h has been computed in [Y, (4.d"), p. 519] as: 



' \(n\) 2 ^k=o ( n -k)\-(n+k)\ J 



In particular the second moment of Ah is thus finite, as was to be shown. What regards 
the even order moments of the Brownian motion B at time h, they are computed as: 

E[Bjr] = ^-T(n + 1/2), 

for any non-negative integer n. The above series thus converges using the ratio test. This 
completes the proof of the normalized price being an entire function in the paramater v. 

The proof of the six-term sum of §6 Theorem being a meromorphic function in v is based 
on the Hermite functions being entire also in their degree. Indeed, from [L, (10.2.8), p. 285] 
one has for any complex numbers li and z the representation: 

2"-r(l/2) / n 1 2 N 2"-r(-l/2) 3 A 

^ ) = r((i-^)/2)'H-2'2'0 + "- r(-^/2) '2'*> 

Herein the reciprocal of the gamma function is entire by construction, and the confluent 
hypergeometric function $ is entire in its first and third variable. 

There is a localization principle for proving entire a function on the complex plane. Indeed, 
by definition one has to show analyticity at any fixed complex number. For this, it is 
sufficient to restrict the function to any relatively compact or compact neighborhood of 
this fixed complex number and prove the function thus obtained analytic. 

What regards the hyperbolic terms, the claim is that as a function of v they can be 
extended as analytic functions from non-negative v to the whole complex plane. For this 
question it is sufficient to consider the function A given by: 

My)= J~ #-(„ + 4) ( K - COSh ^ i9) )E b{u) (h)(y+ie)dy, 
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where b(u) equals ±u or ±(i/+4) and recalling 



Ee(h)(w) = e W? Erfc 




for any complex number w. Extension of this function is essentially by reduction to the case 
where the real part of the degree — (z/+4) of the Hermite function factor in the integrand 
of A(v) is negative, or equivalently, the real part of v is bigger than minus four. 

To fix ideas first consider the case where v is such that the degree of the Hermite function 
factor in the integrand of A{y) is a non-negative integer. This Hermite function is then 
the corresponding Hermite polynomial. In the integrand of A{y) the absolute values of 
the Hermite function factor and the exponential function factor so have linear exponential 
order in the variable y. The decay to zero of the absolute value of the function E b ^ : 
however, is of square exponential order in y. It thus dominates the asymptotic behaviour 
with y to infinity. The absolute value of the integrand of A(u) is so majorized by an 
integrable function, and A can be extended to the above values of v. 

For the general case of the reduction, fix any v of non-positive real part that is not an inte- 
ger. Apply the above localization principle and let v belong to a sufficiently small compact 
neigborhood U in the half-plane {Re (z) < 0}. Shrinking U if necessary assume that for 
any v in U the degree of the Hermite function in the integrand of A{y) is not an integer. 
Using the recursion rule for Hermite functions of §2 express the Hermite function factor of 
A(y) in terms of weighted Hermite functions of negative degrees. Further shrinking U if 
necessary, assume that the so obtained relation represents B on U. Herein the weighting 
factor for the respective Hermite functions are given as powers ofz= -(2g)- 1 / 2 cosh(y+^) 
times polynomials in v. The absolute values of the polynomials in v can be majoriozed 
uniformly on U. The problem thus reduces to majorize by an integrable function on the 
positive real line in the variable y finitely many functions on U times the positive real line 
sending v and y to: 



where k, £ range over finitely many non-negative integers and k is such that is + 4 + k is 
positive. The leading terms of the asymptotic expansion of §2 for any Hermite function 
Hp(z) with degree \i any complex number with negative real part has order z M . Asymp- 
totically with y to infinity, the Hermite function with the smallest positive number u+4+k 
thus dominates the other Hermite function factors in the above functions. With v ranging 
over a compact set, there is a minimal such degree on U. Similarly, there are such ma- 
jorizing choices £* for the factors cosh i (y +i9) and v* for the absolute values of the factors 
Eb(v)(h)(y+id). A three-factor-majorizing function on the positive real line thus results 
whose asymptotic behaviour with y to infinity is governed by the square-exponential decay 
to zero of the corresponding factor E^ u *^(h) and which is integrable. 

If the real part of v is non-negative, the above argument holds in a simplified form. The 
upshot so is that any complex number not an integer has a sufficiently small compact 
neighborhood such that the absolute value of the integrand of A{y) on U times the posi- 
tive real line can be majorized by an integrable function on the positive real line. Herein, 



cosh(y+i9) 



■ cosh\y+iO) ■ E b(iy) (h)(y+i9) , 
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compact neigborhoods can be replaced by relatively compact neigborhoods mutatis mu- 
tandis. Using Lebesgue Dominated Convergence, A can thus be extended as a continuous 
function to the whole complex plane with the integers less than or equal to minus four 
deleted. 

The idea for showing A analytic as a function of v on the complex plane with the integers 
less than or equal to minus four deleted is as follows. Show that differentiation of A is by 
differentation under the integral sign and use that its integrand is entire as function of v. 
For this again first localize to v in any sufficiently small compact neighborhood containing 
no integers less than or equal to minus four. The aim is then to majorize the absolute 
value of the derivative with respect to v of the integrand of A by an integrable function 
independent of v as above. The above argument for getting such a majorizing function 
is based on a comparison of decay rates. The integrand of A has one factor which on 
the positive real line decays to zero of square exponential order whereas the other factors 
explode of at most linear exponential order. This situation is preserved on differentiation 
with respect to the parameter v. In particular, differentiating with respect to the degree 
the asymptotic expansion for Hermite functions on the right half-plane gives an asymptotic 
expansion for this function's partial derivative with respect to the degree. 

At this stage it remains to extend A analytically to the intgers less than or equal to minus 
four. However, A remains bounded in any punctured compact neigbourhood of such an 
integer. Thus A can be extended to an entire function, completing the proof of §5 Theorem. 

Proving the trigonometric terms entire is a direct application of the above localization 
principle to their following factors: 



where h{y) equals v or u+2. In fact, restrict v to belong to any suffiently small compact 
neighborhood U of any point fixed in the complex plane. The above integrand is analytic 
as a function in v and smooth as a function in <fi. Thus the absolute values of its derivatives 
with respect to v of any order are bounded on the product of U and the closed interval 
between zero and 9. In this local situation, using a standard consequence of Lebesgue 
Dominated Convergence, differentation of B with respect to v is thus by partial differen- 
tiation with respect to the parameter v under the integral sign. This proves B entire as a 
function in v. 



16. Trigonometric terms for unit radii: At this point it remains to establish the 
exact form of the trigonometric terms Ctrig,0,6(i/) (0) in §6 Remark for all u, thus extending 
§13 Lemma. This is essentially based on the above integral B being entire. Indeed, B 
being entire allows to first extend the complex identity underlying §6 Remark to all v. 
Thereupon taking real parts then completes the extension argument. 
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